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(57) Zusammenfassung: Die Erfindung betriffl ein Compu- 
tertomographie-Verfahren, bei dem ein Untersuchungsbe- 
reich entlang einer kreisfOrmigen Trajektorie von einem fa- 
cherformigen StrahlenbUndel durchstrahit wird. Im Unter- 
suchungsbereich koharent gestreute Strahlung wird durch 
eine Detektoreinheit gemessen, wobei aus diesen Mess- 
werten der raumliche Verlauf der Streust3ri<e in dem Unter- 
suchungsbereich rekonstruiert wird. Die Rekonstruktion er- 
folgt mittels einer RQckprojekticn in einem Volumen, das 
durch zwei linear unabhangige Vektoren der Rotationsebe- 
ne und einem Wellenvektorubertrag aufgespannt ist. 
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Beschreibung 

[0001] Die Erfindung betrifft ein Computertomogra- 
phie-Verfahren, bei dem ein Untersucfiungsbereich 
entlang einer I<reisf6rmigen Trajel<torie von einem fa- 
cherformigen StraPilenbundel durchstrahit wird und 
eine Detel<toreinheit Primarstrahlung und im Unter- 
suchungsbereich koharent gestreute Strahlen er- 
fasst. Die Erfindung bezielit sicli aulierdem auf einen 
Computertomographen zur DurchfUhrung dieses 
Verfahrens sowie auf ein Computerprogramm zur 
Steuerung des Computertomographen. 
[0002] In bel<annten Verfaliren der eingangs ge- 
nannten Art warden im Untersuchungsbereich koha- 
rent gestreute Strahlen durch eine Detektoreinheit 
gemessen, wobei aus diesen Messwerten der raum- 
liche Verlauf der Streustarke in dem Untersuchungs- 
bereich rekonstruiert wird. Die Rekonstruktion erfolgt 
dabei meist mil iterativen Verfahren basierend auf al- 
gebraischen Rekonstruktionstechniken (ART), bei 
denen eine Rekonstruktion des gesamten durch- 
strahlten Untersuchungsbereichs notwendig ist. Eine 
Rekonstruktion eines Teils des Untersuchungsbe- 
reichs ist nicht moglich. Aufierdem erfordern die ite- 
rativen Verfahren einen hohen Rechenaufwand und 
fuhren somit zu relativ langen Rekonstruktionszeiten. 
Alternative Verfahren sind zweidimensionale, auf 
Ruckprojektion basierende Rekonstruktionstechni- 
ken, bei denen zwischen dem Untersuchungsbereich 
und der Detektoreinheit so grof^ Abstdnde bestehen 
bzw. angenommen werden, dass eine eindeutige Be- 
ziehung zwischen dem AuftrefFpunkt eines Streu- 
strahls auf die Detektoreinheit und dem Streuwinkel 
vorhanden ist. Diese Geometrie kann jedoch in medi- 
zinischen Anwendungen nicht realisiert werden, so 
dass die zweidimensionale Ruckprojektion nur unge- 
naue Ergebnisse liefert. 

[0003] Es ist daher Aufgabe der vorliegenden Erfin- 
dung, ein Verfahren, einen Computertomographen 
sowie ein Computerprogramm anzugeben, mit denen 
eine schnellere Rekonstruktion der Streustarkever- 
teilung im Untersuchungsbereich unter Berucksichti- 
gung der tatsdchlichen Geometrie in medizinischen 
Anwendungen mOglich ist. 

[0004] Diese Aufgabe wird in verfahrensm3(liger 
Hinsicht erfindungsgemflli durch ein Computertomo- 
graphie- Verfahren gelost mit den Schritten 

- Erzeugen eines einen Untersuchungsbereich 
bzw. ein darin befindliches Objekt durchsetzen- 
den Strahlenfachers mit einer Strahlenquelle und 
einer Blenden-Anordnung, die sich zwischen dem 
Untersuchungsbereich und der Strahlenquelle be- 
findet, 

- Erzeugen von Relativbewegungen zwischen der 
Strahlenquelle einerseits und dem Untersu- 
chungsbereich bzw. dem Objekt andererseits, die 
eine Rotation urn eine Rotationsachse umfassen, 

- Akquisition von Messwerten, die abhangig von 
der Intensitat der Strahlung sind, mit einer Detek- 
toreinheit, die Primarstrahlung aus dem Strahlen- 



facher und in dem Untersuchungsbereich bzw. an 
dem Objekt kohdrent gestreute Strahlung erfasst, 
wahrend der Relativbewegungen, 
- Rekonstruktion eines CT-Bildes des Untersu- 
chungsbereichs aus den Messwerten, wobei eine 
Ruckprojektion in einem Volumen, das durch zwei 
linear unabhangige Vektoren der Rotationsebene 
und einen WellenvektorUbertrag aufgespannt ist, 
durchgefuhrt wird. 

[0005] Die Streustdrke hSngt nicht allein vom Mate- 
rial, sondern auch von dem Streuwinkel und der Wel- 
lenlange der Strahlung ab. Eine Ruckprojektion in ei- 
nem Volumen, das durch zwei linear unabhangige 
Vektoren der Rotationsebene und einen Wellenvek- 
torubertrag aufgespannt ist, hat dagegen den Vorteil, 
dass die so parametrisierte Streustarke nur noch vom 
Streumaterial abhangig ist, da namlich der Wellen- 
vektorUbertrag bekanntlich proportional zum Produkt 
aus der inversen Welleniange und dem Sinus des 
halben Streuwinkels ist. Der Streuwinkel ist hierbei 
der Winkel, den der Verlauf des gestreuten Strahls 
mit dem Verlauf einschliefit, den der Strahl ohne den 
Streuprozess genommen hatte. 
[0006] In dem oben genannten Volumen weisen die 
Streustrahlen eine gekrummte Form auf. Die Beruck- 
sichtigung dieser gekrummten Form der Streustrah- 
len bei der Ruckprojektion gemafi Anspruch 2 fuhrt 
zu einer verbesserten Qualitat der rekonstruierten 
Stre usta rkeverteil u ng . 

[0007] In dem bevorzugten Rekonstruktionsverfah- 
ren gemaH Anspruch 3 werden die Messwerte vor 
der Ruckprojektion gewichtet. Dabei berUcksichtigt 
die Wichtung die abnehmende effektive Detektorfla- 
che mit zunehmenden Streuwinkel und die abneh- 
mende Strahldichte mit zunehmenden Abstand des 
Auftreffpunkts des gestreuten Strahles auf der Detek- 
toreinheit vom Streuzentrum. Die effektive Detektor- 
flache ist die in eine Ebene senkrecht zum Verlauf 
des detektieren Strahles projizierte Detektorflache. 
Die Wichtung im Anspruch 4 berUcksichtigt die ab- 
nehmende Strahldichte mit zunehmenden Abstand 
der Strahlenquelle vom Streuzentrum. Die BerUck- 
sichtigung dieser EfFekte fUhrt zu einer verbesserten 
Qualitat der rekonstruierten Streustarkeverteilung. 
[0008] Anspruch 5 beschreibt ein bevorzugtes Re- 
konstruktionsverfahren mit einem im Vergleich zu an- 
deren Methoden geringeren Rechenaufwand, das zu 
einer guten Bildqualitat fuhrt. Ein Computertomo- 
graph zur Durchfuhrung des erfindungsgemSlien 
Verfahrens ist in Anspruch 6 beschrieben. 
[0009] Anspruch 7 definiert ein Computerprogramm 
zur Steuerung eines Computertomographen nach 
Anspruch 6. 

[0010] Die Erfindung wird nachfolgend anhand ei- 
nes Ausfuhrungsbeispiels eriautert, wobei auf eine 
Zeichnung Bezug genommen wird, in der 
[0011] Fig. 1 einen Computertomographen zeigt, 
mit dem das erfindungsgemaiie Verfahren ausfUhr- 
bar ist, 
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[0012] Fig. 2 eine schematische Querschnittsan- 
sicht des Computertomographen nach Fig. 1 in Rich- 
tung der Rotationsachse zeigt, 
[0013] Fig. 3 ein Ablaufdiagramm einer Ausfuh- 
rungsform des erfindungsgemaBen Verfahrens zeigt, 
[0014] Fig. 4 eine Draufsiclit auf eine Spalte von 
Detektorelementen zeigt, 

[0015] Fig. 5 eine schematische Darstellung der 
Anordnung virtueller Strahlenquellen zeigt, 
[0016] Fig. 6 eine Schnittfiache durch die Stralilen 
der virtuellen Strahlenquellen zeigt, 
[0017] Fig. 7 eine schematische Schnittansicht des 
durchstrahlten Untersuchungsbereichs und des De- 
tektors in Richtung der Rotationsebene zeigt, und 
[0018] Fig. 8 die Abhangigkeit des Betrags des 
Wellenvektorubertrags von dem Abstand eines 
Streuzentrums im Untersuchungsbereich von dem 
Fudpunkt des Detektors (Auftreffpunkt des Primar- 
strahls auf die Detektoreinheit) zeigt. 
[0019] Der in Fig. 1 dargestellte Computertonno- 
graph umfasst eine Gantry 1, die um eine Rotations- 
achse 14 rotieren kann. Dazu wird die Gantry 1 von 
einem Motor 2 mit einer vorzugsweise konstanten, 
aber einstellbaren Winkelgeschwindigkeit angetrie- 
ben. An der Gantry 1 ist eine Strahlenquelle S befes- 
tigt, beispielsweise ein Rontgenstrahler. Eine Blen- 
denanordnung 31 legt einen zur Untersuchung be- 
nutzten Strahlenfacher 41 fest, der in Fig. 1 mit 
durchgezogenen Linien dargestellt ist. Der Strahlen- 
facher 41 veriauft senkrecht zur Rotationsachse 14 
und hat in deren Richtung geringe Abmessungen, 
beispielsweise 1 mm. Zwischen der Blendenanord- 
nung 31 und der Strahlenquelle S kann eine zweite 
Blendenanordnung 32 angeordnet sein, die aus der 
von der Strahlenquelle S erzeugten Strahlung ein ke- 
gelformiges Strahlenbundel 42 ausblendet. Das ke- 
gelformige Strahlenbundel 42, das ohne die Blenden- 
anordnung 31 entstehen wurde, ist durch gestrichelte 
Linien angedeutet. 

[0020] Der Strahlenfacher 41 durchdringt einen zy- 
linderformigen Untersuchungsbereich 13, in dem 
sich ein Objekt, z. B. ein Patient auf einem Patienten- 
lagerungstisch (beides nicht dargestellt) oder aber 
auch ein technisches Objekt befinden kann. Nach 
dem Durchsetzen des Untersuchungsbereichs 13 
trifft der Strahlenfacher 41 auf eine an der Gantry 1 
befestigte Detektoreinheit 16 mit einer Messflache 
17, die eine Vielzahl von matrixformig angeordneten 
Detektorelementen umfasst. Die Detektorelemente 
sind in Zeilen und Spalten angeordnet. Die Detektor- 
spalten verlaufen parallel zur Rotationsachse 14. Die 
Detektorzeilen befinden sich in zur Rotationsachse 
senkrechten Ebenen, vorzugsweise auf einem Kreis- 
bogen um die Strahlenquelle S. Sie kOnnen aber 
auch anders geformt sein, z. B. einen Kreisbogen um 
die Rotationsachse 14 beschreiben oder geradlinig 
sein. Die Detektorzeilen enthalten in der Regel we- 
sentlich mehr Detektorelemente (z.B. 1000) als die 
Detektorspalten (z.B. 16). 

[0021] Der Strahlenfacher 41, der Untersuchungs- 



bereich 13 und die Detektoreinheit 16 sind aneinan- 
der angepasst. In einer zur Rotationsachse 14 senk- 
rechten Ebene sind die Abmessungen des Strahlen- 
fachers 41 so gewahit, dass der Untersuchungsbe- 
reich 13 vollstandig durchstra hit wird, und die Lange 
der Zeilen der Detektoreinheit 16 ist gerade so be- 
messen, dass der Strahlenfacher 41 vollstandig er- 
fasst werden kann. Der Strahlenfacher trifft auf den 
bzw. auf die mittlere(n) Detektorzeile(n). 
[0022] Wenn es sich um ein technisches Objekt 
handelt und nicht um einen Patienten, kann das Ob- 
jekt bei einer Untersuchung gedreht werden, wah- 
rend die Strahlenquelle S und die Detektoreinheit 16 
stillstehen. Der Untersuchungsbereich 13 bzw. das 
Objekt Oder der Patientenlagerungstisch kann mittels 
eines Motors 5 parallel zur Rotationsachse 14 ver- 
schoben werden. Dazu aquivalent konnte aber auch 
die Gantry in diese Richtung verschoben werden. 
[0023] Wenn die Motoren 2 und 5 gleichzeitig lau- 
fen, beschreiben die Strahlenquelle S und die Detek- 
toreinheit 16 eine helixfOrmige Trajektorie relativ zum 
Untersuchungsbereich 13. Wenn hingegen der Motor 
5 fur den Vorschub in Richtung der Rotationsachse 
14 stillsteht und der Motor 2 die Gantry rotieren lasst, 
ergibt sich eine kreisformige Trajektorie fur die Strah- 
lenquelle S und die Detektoreinheit 16 relativ zum 
Untersuchungsbereich 13. Im folgenden wird nur die 
kreisformige Trajektorie betrachtet. 
[0024] Fig. 2 zeigt, dass sich zwischen dem Unter- 
suchungsbereich 13 und der Detektoreinheit 16 eine 
Kollimator-Anordnung 6, die eine Vielzahl ebener La- 
mellen 60 umfasst, befindet. Die Lamellen 60 beste- 
hen aus einem die Rdntgenstrahlung stark absorbie- 
renden Material und liegen in Ebenen, die parallel zur 
Rotationsachse 14 verlaufen und sich im Fokus der 
Strahlenquelle S schneiden. Ihr Abstand kann bei- 
spielsweise 1 cm betragen, und jede Lamelle 60 kann 
in der Zeichenebene eine Abmessung von z.B. 20 cm 
haben. Durch die Kollimator-Anordnung 6 wird also 
der Strahlenfacher 41 in eine Anzahl von einander 
benachbarten Abschnitten unterteilt, so dass eine 
Spalte von Detektorelementen im wesentlichen nur 
von PrimSr- bzw. Streustrahlen aus einem Abschnitt 
getroffen wird. 

[0025] Die von der Detektoreinheit 16 akquirierten 
Messwerte werden einem Bildverarbeitungsrechner 
lOzugefuhrt, der mit der Detektoreinheit 16 z. B. uber 
eine kontaktlos arbeitende Datenubertragung (nicht 
dargestellt) verbunden ist. Der Bildverarbeitungs- 
rechner 10 rekonstruiert die Streustarkeverteilung im 
Untersuchungsbereich 13 und gibt sie, beispielswei- 
se auf einem Monitor 11, wieder. Die beiden Motoren 
2 und 5, der Bildverarbeitungsrechner 10, die Strah- 
lenquelle S und der Transfer der Messwerte von der 
Detektoreinheit 16 zum Bildverarbeitungsrechner 10 
werden von einer Kontrolleinheit 7 gesteuert. 
[0026] In anderen Ausfuhrungsformen konnen die 
akquirierten Messwerte zur Rekonstruktion zunachst 
einem oder mehreren Rekonstruktionsrechnern zu- 
gefuhrt werden, die die rekonstruierten Daten z.B. 
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Uber ein Glasfaserkabel an den Bildverarbeitungs- 
rechner weiteiieiten. 

[0027] Fig. 3 zeigt den Ablauf einer Ausfuhrungs- 
form eines Mess- und Rekonstruktionsverfahrens, 
das mit dem Computertomographen nach Fig. 1 
durchgefuhrt werden kann. 

[0028] Nach der Initialisierung im Schritt 101 rotiert 
die Gantry mit einer konstanten Winkelgescliwindig- 
keit. Im Schritt 103 wird die Strahlung der Strahlen- 
quelle S eingeschaltet, so dass die Detektoreinheit 
16 die Prinndr- und die Streustrahlung aus einer Viel- 
zahl von Winkelpositionen erfassen kann. Das oder 
die Detektorelemente in der Mitte jeder Detektorspal- 
te erfassen im wesentlichen die Primarstrahlung, 
wahrend die Streustrahlung (Sekundarstrahlung) von 
den in einer Spalte v^eiter auBen liegenden Detektor- 
elemente n erfasst wird. 

[0029] Dies ist in Fig. 4 schematisch angedeutet, 
die eine Draufsicht auf eine Spalte von Detektorele- 
menten darstellt. Die Detektorelemente 161, die die 
Streustrahlung erfassen, sind einfach schraffiert dar- 
gestellt, wahrend das Detektorelement 160 im Zen- 
trum, das die Primarstrahlung erfasst, durch ein 
Kreuz markiert ist. Bei anderen Strahlenquellen, ins- 
besondere bei Strahlenquellen mit einem grolJeren 
Fokus, konnten auch mehr als ein Detektorelement 
die Primarstrahlung erfassen. Beiderseits dieses 
zentralen Detektorelements befinden sich auf Grund 
der endlichen Abmessungen des Fokus der Strahlen- 
quelle Detektorelemente, die von Streustrahlung, 
aber auch von einer (reduzierten) PrimSirstrahlung 
getrofFen werden (weif^e Elemente benachbart zu 
160). In dieser Ausfuhrungsform werden als Streu- 
strahlen nur die Strahlen angesehen, die von den in 
der Zeichnung schraffiert dargestellten Detektorele- 
menten gemessen werden. 

[0030] Die Streustarke ist unter anderem abhangig 
von der Energie des gestreuten Rontgenquants. Da- 
her muss die Energie des gestreuten Rontgenquants 
entweder gemessen werden, was voraussetzt, dass 
die Detektorelemente Energie aufldsend messen 
konnen, oder es muss Rdntgenstrahlung mit Quante- 
nenergien aus einem mdglichst kleinen Bereich (im 
Idealfall monochromatische Rdntgenstrahlung) ver- 
wendet werden. Um die Energiedifferenz der ROnt- 
genquanten relativ zur ihrer Energie mdglichst zu mi- 
nimieren, gibt es verschiedene Moglichkeiten: 

- Der Gebrauch geeigneter Filtermaterialien, z.B. 
Kupfer, im Primarstrahl. Dadurch wird die von ei- 
nem Rontgenstrahler erzeugte weiche Rontgen- 
strahlung, d.h. Rontgenstrahlung mit niedriger 
Quantenenergie, weitgehend unterdruckt. 

- Zusatzlich kann die Spannung einer Rdntgen- 
rOhre bezUglich des gewdhlten Filters optimiert 
werden. 

- Schliel^lich ist es moglich, die sogenannte 
„balnced filter"-Technik anzuwenden. Dabei wer- 
den die Daten zweimal akquiriert, wobei sich im 
Strahlengang jeweils Filter mit leicht voneinander 
abweichender Ordnungszahl befinden, deren 



K-Kante zur Filterung herangezogen wird. An- 
schliedend wird das Diffierenzsignal aus den bei- 
den Messungen extrahiert. 

[0031] Im Schritt 105 werden die Messwerte der 
Streustrahlen normiert. Dabei werden die Messwerte 
jeder Strahlenquellenposition der gestreuten Strah- 
len durch die Messwerte derjenigen Primarstrahlen 
dividiert, die die gestreuten Strahlen verursacht ha- 
ben. 

[0032] Im Schritt 107 erfolgt ein Rebinning der 
Messwerte. Dabei wird jedem Messwert eine Linie 
von dem Detektorelement, an dem der Messwert de- 
tektiert worden ist, zur Strahlenquellenposition zuge- 
ordnet. Es wird also angenommen, dass Strahlen des 
fiktiven, kegelformigen Strahlenbundels 42 die Mess- 
werte verursacht haben, ohne dass die Strahlen ge- 
streut worden waren. Durch das Rebinning werden 
die Messwerte nun so umsortiert, als wdren sie mit ei- 
ner anderen Strahlenquelle (einer kreisbogenfOrmi- 
gen Strahlenquelle, die jeweils zueinander parallele 
Strahlenfacher emittieren kann) und mit einem ande- 
ren Detektor (einem ebenen, rechteckigen virtuellen 
Detektor) gemessen worden. 

[0033] Dies wird anhand von Fig. 5 naher eriautert. 
Mit 17 ist dabei die kreisformige Trajektorie bezeich- 
net, von der aus die Strahlenquelle den Untersu- 
chungsbereich durchstrahlt. Mit 413 ist ein facherfor- 
miges Strahlenbundel bezeichnet, das von der Strah- 
lenquellenposition So ausgeht und dessen Strahlen 
in einer die Rotationsachse 14 enthaltene Ebene ver- 
laufen. 

[0034] Man kann sich das kegelfdrmige Strahlen- 
bundel, das von der Strahlenquellenposition S^ emit- 
tiert wird, aus einer Vielzahl von ebenen Strahlenfa- 
chern zusammengesetzt denken, die sich in zur Ro- 
tationsachse 14 parallelen Ebenen befinden, die sich 
in der Strahlenquellenposition schneiden. Fig. 5 zeigt 
von diesen Strahlenfachern nureinen einzigen, nSm- 
lich den Strahlenfacher 413. 

[0035] Auderdem sind in Fig. 5 noch weitere Strah- 
lenfacher 411, 412 und 414, 415 dargestellt, die par- 
allel zu dem Strahlenfacher 413 und in zueinander 
und zur Rotationsachse 14 parallelen Ebenen liegen. 
Die zugehorigen Strahlenquellenpositionen S_2, S_, 
und S^, Sj werden von der Strahlenquelle S einge- 
nommen, bevor bzw. nachdem sie die Strahlenquel- 
lenposition So erreicht hat. Alle Strahlen in den Strah- 
lenfachern 411 bis 41 5 haben denselben Projektions- 
winkel. Der Projektionswinkel ist der Winkel, den die 
Ebenen der Strahlenfacher mit einer zur Rotations- 
achse 14 parallelen Referenzebene einschliel^en. 
[0036] Die Strahlenfacher 411 bis 41 5 definieren ein 
Strahlenbundel 410 mit einer zeltartigen Form. In 
Fig. 5 und 6 ist der Schnittbereich 420 dargestellt, 
der sich ergibt, wenn das Strahlenbundel 410 von ei- 
ner die Rotationsachse 14 enthaltenden und zu den 
Ebenen der Strahlenfacher 411 bis 415 senkrechten 
Ebene geschnitten wird. Die beiden die Rotations- 
achse schneidenden Rander des Schnittbereichs 
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420 sind gewOlbt. Diese W6lbung ist darauf zurQck- 
zufUhren, dass die Strahlenquellenpositionen in der 
Mitte (z.B. Sq) weiter von der Schnittebene entfernt 
sind als die am Rande (z.B. oder S_2) und dass die 
Faclier alle den gleiclien Offnungswinl<el aufweisen. 
Fur jede Gruppe von Strahlenfachern wird daher in 
dem ebenen Schnittbereich 420 ein rechteckiger vir- 
tueller Detektor 170 definiert, dessen Rander 171 
bzw. 1 72 durcli die Abmessungen der aufieren Strali- 
Ienfdcher411 bzw. 415 in dem ebenen Schnittbereich 
gegeben sind. 

[0037] In Fig. 6 sind auBerdem - durch runde Punk- 
te markiert - die Durchstolipunkte einiger in den 
Strahlenfachern 411...415 enthaltener Strahlen durch 
diesen virtuellen Detektor dargestellt. Schliefilich 
sind durch Kreuze die Stutzpunkte eines regelmafli- 
gen kartesischen Gitters angedeutet. Die Durchstoli- 
punkte und die Stutzpunkte fallen durchweg nicht zu- 
sammen. Aus den Messwerten fOr die DurchstoH- 
punkte mUssen daher die Messwerte an den Squidis- 
tanten StUtzpunkten innerhalb des virtuellen Detek- 
tors 170 ermittelt werden. 

[0038] Im einzelnen ist dieses Rebinning in der DE 
198 45 133 A1 beschrieben, worauf hiermit Bezug 
genommen wird. 

[0039] Im Schritt 109 wird auf die aus dem Rebin- 
ning hervorgehenden Messwerte der Streustrahlen 
eine eindimensionale Filterung mit einem rampenfor- 
mig mit der Raumfrequenz ansteigenden Ubertra- 
gungsfaktor angewandt. Dazu werden jeweils aufein- 
ander folgende Werte in einer Richtung parallel zur 
Rotationsebene, also entlang einer Zeile des virtuel- 
len Detektors, herangezogen. Diese Filterung wird 
entlang jeder Zeile des virtuellen Detektors fur alle 
Projektionswinkel durchgefuhrt. 
[0040] In anderen Ausfuhrungsformen konnte auf 
das Rebinning verzichtet werden. Dann ist bekannt- 
lich, da die Detektoreinheit z.B. bogenformig um die 
Strahlenquelle oder um die Rotationsachse ge- 
krummt ist, die Filterung zu modifizieren. 
[0041] Die gefilterten Messwerte werden auschlie- 
l^end zur Rekonstruktion der Streustarkeverteilung 
im Untersuchungsbereich durch eine RUckprojektion 
herangezogen. 

[0042] Die RQckprojektion erfolgt dabei in einem Vo- 
lumen, das durch die Vektoren x, y und q aufge- 
spannt ist, wobei die Einheitsvektoren x und y in der 
Rotationsebene liegen und horizontal bzw. vertikal 
orientiert sind und der Wellenvektorubertrag q paral- 
lel zur Rotationsachse orientiert ist. In anderen Aus- 
ftihrungsbeispielen konnen statt der Vektoren x und 
y zwei andere, linear unabhangige Vektoren der Ro- 
tationsebene verwendet werden. Der Betrag des 
WellenvektorUbertrags q ist, wie oben bereits er- 
wahnt, proportional zum Produkt der inversen Wel- 
leniange A der gestreuten ROntgenquanten und des 
Sinus des halben Streuwinkels 0 

q = (1/A)sin(0/2). (1) 



[0043] Der Streuwinkel O lasst sich anhand der in 
Fig, 7 gezeigten Anordnung des vom FScherstrahl 
41 durchstrahlten Untersuchungsbereichs 15 und der 
Detektoreinheit 16 bestimmen. In dieser Figur wird, 
im Gegensatz zu bekannten zweidimensionalen, auf 
RQckprojektion basierenden Verfahren, die reale Ge- 
ometrie in medizinischen Anwendungen berucksich- 
tigt. Dies fuhrt fur ein bzw. jedes Detektorelement Dj 
zu unterschiedlichen Streuwinkein 0, die gemali fol- 
gender Gleichung berechnet werden: 

0 = arctan(a/d ). (2) 

[0044] Dabei bezeichnen d den Abstand eines 
Streuzentrums S, und a den Abstand des Detektore- 
lements D^ von dem Fulipunkt 12 des Detektors. 
[0045] Das Detektorelement D^ erfasst Strahlen, die 
im vom Facherstrahl 41 durchstrahlten Untersu- 
chungsbereich 15 unter den Winkein 0^ < © < ©2 ge- 
streut werden. 

[0046] Aus den beiden oben genannten Gleichun- 
gen ergibt sich fur kleine Winkel © 

q = a/(2dA). (3) 

[0047] Der Betrag des Wellenvektorubertrags q in 
Abhangigkeit vom Abstand d eines Streuzentrums 
vom Fufipunkt des Detektors hat somit einen hyper- 
belartigen und damit nichtlinearen Verlauf (darge- 
stellt in Fig. 8). Daraus folgt, dass der ursprunglich 
gerade Verlauf der Streustrahlen im (x,y,q)-Raum ge- 
krUmmt ist. Die RQckprojektion erfolgt also entlang 
hyperbelartig gekriimmter Strahlen. 
[0048] Im Schritt 111 wird ein Voxel V(x,y,q) inner- 
halb eines vorgebbaren (x,y) - Bereiches (field of 
view - FOV) und innerhalb eines sich aus der Geo- 
metrie des Computertomographen ergebenden Wer- 
te be re ichs des Wellenvektorubertrags q bestimmt. 
[0049] Im Schritt 113 werden die gefilterten Werte 
mit einem Gewichtungsfaktor multipliziert, der dem 
Kehrwert des Cosinus des Streuwinkels entspricht. 
Dadurch wird die abnehmende effektive Detektorfla- 
che mit zunehmendem Streuwinkel berOcksichtigt. 
Wenn der genannte Winkel klein ist, ist der Cosinus 
des Winkels praktisch stets 1 , so dass diese Wich- 
tung entfallen kann. Aufierdem wird die abnehmende 
Strahldichte mit zunehmendem Abstand des Streu- 
zentrums, d.h. des Voxels V(x,y,q), vom Auftreffpunkt 
des Streustrahls auf die Messflache berucksichtigt, 
indem alle Messwerte fur jede Strahlenquellenpositi- 
on mit einem Gewichtungsfaktor multipliziert werden, 
der dem Quadrat des Abstandes zwischen dem 
Streuzentrum, an dem der zum Messwert gehSrende 
Strahl gestreut wurde, und dem Auftreffpunkt des 
Streustrahls auf die Messflache entspricht. 
[0050] Wenn in anderen Ausfuhrungsformen auf 
das Rebinning verzichtet wird, so ist eine zusatzliche 
Multiplikation der gefilterten Messwerte mit einem 
Gewichtungsfaktor erforderlich, der dem Kehrwert 
des Quadrats des Abstandes zwischen der Strahlen- 
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quellenposition und dem Streuzentmm entspricht, an 
dem der detektierte Strahl gestreut wurde. 
[0051 ] Bei der Ruckprojektion im Schritt 115 werden 
nun alle gekrunnmten Strahlen herangezogen, die 
den Voxel V(x,y,q) passieren. Passiert kein Strahl ei- 
ner Strahlenquellenposition exakt die Mitte des Vo- 
xels, dann muss der zugehorige Wert durch Interpo- 
lation der Messwerte von benachbarten Strahlen er- 
mittelt werden. Der Messwert, der dem den Voxel 
passierenden Strahl zugeordnet werden kann, bzw. 
der durch Interpolation gewonnene Messwert wird 
auf den Voxel V(x,y,q) akkumuliert. Nachdem auf die- 
se Weise fur alle Strahlenquellenpositionen die Bei- 
trage zu dem betreffenden Voxel akkumuliert worden 
sind, wird im Schritt 117 gepruft, ob alle Voxel in dem 
zu rekonstruierenden (x,y,q) - Bereich durchlaufen 
sind. 1st dies nicht der Fall, verzweigt das Ablaufdia- 
gramm zum Schritt 111. Ansonsten ist die Streustar- 
keverteilung fUr alle Voxel in dem FOV ermittelt wor- 
den, und das Rekonstruktionsverfahren ist beendet 
(Schritt 119). 

[0052] In anderen Ausfuhrungsformen kann die 
Ruckprojektion in dem durch zwei linear unabhangi- 
ge Vektoren der Rotationsebene und durch den Wel- 
lenvektorubertrag aufgespannten Volumen als Nahe- 
rung entlang gerader Strahlen durchgefUhrt werden. 

Patentanspruche 

1. Computertomographie-Verfahren mit den 
Schritten: 

a) Erzeugen eines einen Untersuchungsbereich (13) 
bzw. ein darin befindliches Objekt durchsetzenden 
Strahlenfachers (41) mit einer Strahlenquelle (S) und 
einer Blenden-Anordnung (31), die sich zwischen 
dem Untersuchungsbereich (13) und der Strahlen- 
quelle (S) befindet, 

b) Erzeugen von Relativbewegungen zwischen der 
Strahlenquelle (S) einerseits und dem Untersu- 
chungsbereich (13) bzw. dem Objekt andererseits, 
die eine Rotation um eine Rotationsachse (14) um- 
fassen, 

c) Akquisition von Messwerten, die abhSngig von der 
Intensitat der Strahlung sind, mit einer Detektorein- 
heit (16), die Primarstrahlung aus dem Strahlenfa- 
cher (41) und in dem Untersuchungsbereich (13) 
bzw. an dem Objekt koharent gestreute Strahlung er- 
fasst, wahrend der Relativbewegungen, 

d) Rekonstruktion eines CT-Bildes des Untersu- 
chungsbereichs (13) aus den Messwerten, wobei 
eine Ruckprojektion in einem Volumen, das durch 
zwei linear unabhangige Vektoren der Rotationsebe- 
ne und einen WellenvektorUbertrag aufgespannt ist, 
durchgefUhrt wird. 

2. Computertomographie-Verfahren nach An- 
spruch 1, wobei im Rekonstruktionsschritt d) die 
Ruckprojektion entlang von Strahlen gekrummter 
Form durchgefOhrt wird. 



3. Computertomographie-Verfahren nach An- 
spruch 1, wobei im Rekonstruktionsschritt d) die 
Messwerte vor der Ruckprojektion mit einem ersten 
Gewichtungsfaktor, der dem Quadrat des Abstandes 
zwischen dem Streuzentrum, an dem der detektierte 
Strahl gestreut wurde, und dem Auftreffpunkt des 
Streustrahls auf der Detektoreinheit entspricht, und 
mit einem zweiten Gewichtungsfaktor, der dem Kehr- 
wert des Cosinus des Streuwinkels entspricht, multi- 
pliziert werden. 

4. Computertomographie-Verfahren nach An- 
spruch 3, wobei im Rekonstruktionsschritt d) vor der 
Ruckprojektion alle Messwerte fur jede Strahlenquel- 
lenposition mit einem Gewichtungsfaktor, der dem 
Kehrwert des Quadrats des Abstandes zwischen der 
Strahlenquellenposition und dem Streuzentrum, an 
dem der detektierte Strahl gestreut wurde, entspricht, 
multipliziert werden. 

5. Computertomographie-Verfahren nach An- 
spruch 1 , wobei der Rekonstruktionsschritt d) folgen- 

de Schritte umfasst: 

- eindimensionale Filterung der Messwerte in Rich- 
tung parallel zur Rotationsebene, 

- Rebinning der Messwerte zu einer Anzahl von 
Gruppen, wobei jedem von einem Detektorelement 
gemessenen Messwert eine Linie von dem Detektor- 
element zur Strahlenquellenposition zugeordnet wird 
und jede Gruppe mehrere zueinander und zur Rota- 
tionsachse parallele Ebenen umfasst, in denen sich 
je ein Linienfacher (411. ..415) befindet, 

- Rekonstruktion der Verteilung der Streustarke aus 
den Messwerten, wobei eine Ruckprojektion in einem 
Volumen, das durch zwei linear unabhangige Vekto- 
ren der Rotationsebene und einen Wellenvektoruber- 
trag aufgespannt ist, durchgefUhrt wird. 

6. Computertomograph zur DurchfUhrung des 
Verfahrens nach Anspruch 1 mit 

- einer Strahlenquelle (S) und einer Blenden-Anord- 
nung (31), die sich zwischen dem Untersuchungsbe- 
reich (13) und der Strahlenquelle (S) befindet, zum 
Erzeugen eines einen Untersuchungsbereich (13) 
bzw. ein darin befindliches Objekt durchsetzenden 
Strahlenfachers (41), 

- einer mit der Strahlenquelle (S) gekoppelten Detek- 
toreinheit (16) mit einer Messflache (17), 

- einer Antriebsanordnung (2,5), um ein im Untersu- 
chungsbereich (13) enthaltenes Objekt und die 
Strahlenquelle (S) relativ zueinander um eine Rotati- 
onsachse (14) und/oder sich parallel zur Rotations- 
achse (14) verschieben zu lassen, 

- einer Rekonstruktionseinheit (10) zur Rekonstrukti- 
on der Verteilung der Streustarke innerhalb des Un- 
tersuchungsbereichs aus den von der Detektorein- 
heit (16) akquirierten Messwerten, 

- einer Steuereinheit (7) zur Steuerurg der Strahlen- 
quelle (S), der Detektoreinheit (16), der Antriebsan- 
ordnung (2,5) und der Rekonstruktionseinheit (10) 
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entsprechend den Schritten a) bis d) nach Anspruch 
1. 

7. Computerprogramm fur eine Steuereinheit (7) 
zur Steuerung einer Strahlenquelle (S), einer Blen- 
den-Anordnung (31), einer Detektoreinheit (16), einer 
Antriebsanordnung (2,5) und einer Rekonstruktions- 
einheit (10) eines Computertomographen mit einem 
Ablauf nach Anspruch 1 . 

Es folgen 7 Blatt Zeichnungen 
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Fig. 1 




8/14 



DE 102 52 662 A1 2004.05.27 



Fig. 2 
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Fig. 4 
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Fig, 6 
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